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Abstract 

This paper is devoted to investigate the cold plasma wave proper- 
ties outside the event horizon of the Schwarzschild planar analogue. 
The dispersion relations are obtained from the corresponding Fourier 
analyzed equations for non-rotating and rotating, non-magnetized and 
magnetized backgrounds. These dispersion relations provide complex 
wave numbers. The wave numbers are shown in graphs to discuss the 
nature and behavior of waves and the properties of plasma lying in 
the vicinity of the Schwarzschild event horizon. 

Keywords : Schwarzschild planar analogue, cold plasma. 



1 Introduction 

General Relativity is the geometric theory of gravity which deals with four- 
dimensional spacetime. According to this theory, gravity severely modifies 
space and time near the black hole. The Schwarzschild black hole, a minimum 
configuration of the Kerr black hole, exhibits a strong gravitational field 
[Ij. This field pulls the plasma surrounding the black hole event horizon 
towards and then into the black hole. The plasma is drawn equally from all 
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directions and thus form an omni-directional accretion disk. The existence 
of this disk is one of the observational facts used to identify a black hole. 
The motion of plasma inside the accretion disk is governed by the theory of 
general relativistic magnetohydrodynamics (GRMHD). Maxwell's equations, 
Ohm's law, mass, momentum and energy conservation equations constitute 
this theory that are required to investigate various aspects of the interaction 
of relativistic gravity with plasma's magnetic field. 

The theory of general relativity is a sterile subject until it touches the 
real physical world. Only the hard reality of experiments and of astronom- 
ical observations can bring the theory to life. The 3+1 ADM formalism, 
developed to study the quantization of gravitational field by Arnowitt et al. 
[2], has mostly been used in numerical relativity [3]-[Z|. The 3+1 spacetime 
split in the formulation of general relativity is particularly appropriate for 
applications to the black hole theory as described by Thorne et al. |8]-[T0]. 
Wave properties in the Friedmann universe were investigated [II]-[I3] by the 
same formalism. Sakai and Kawata [H] developed a linearized treatment of 
relativistic plasma waves in analogy with the special relativistic formulation. 
Khanna [T^ derived Ohm's law for two component plasma theory of the Kerr 
black hole. Zhang [16]- [T7] formulated stationary symmetric GRMHD theory 
with its applications in Kerr geometry. Buzzi et al. [I8]-[T9] treated the wave 
propagation in radial direction close to the Schwarzschild horizon in general 
relativistic two component plasma. 

In a recent paper, Sharif and Umber [20] have investigated real wave 
numbers and properties of the medium existing in the vicinity of a Schwarz- 
schild black hole using the same split. The cold plasma is considered in 
non-rotating or rotating, non-magnetized or magnetized backgrounds. The 
same authors have also worked out the wave properties in isothermal plasma 
for the Schwarzschild spacetime planar analogue using ADM split [2I]-[22]- It 
is verified that no information can be extracted either from the event horizon 
or its exterior. 

In this paper, we focus our attention to investigate the wave numbers and 
properties of the medium for the Schwarzschild planar analogue spacetime. 
The paper is organized as follows. Section 2 describes the general line element 
and its specified form (the planar analogue) for the Schwarzschild black hole 
in 3+1 formalism. It also contains background assumptions. In sections 
3, 4 and 5, the dispersion relations for non-rotating (non- magnetized and 
magnetized), rotating non- magnetized and rotating magnetized backgrounds 
are investigated respectively. Finally, section 6 provides the outlook of the 



2 



results. 



2 3+1 Spacetime Modeling 

The general line element in 3 + 1 formalism can be expressed as [T7] 

ds"^ = -a^df + -fijidx' + (5'dt){dx^ + (3Ut), 

where a is the lapse function, are the components of the shift vector and 
are the components of the spatial metric. All these quantities are functions 
of coordinates t and x*. A natural observer, associated with this spacetime 
called fiducial observer (FIDO), has four- velocity perpendicular to the hy- 
persurfaces of constant time t. The planar analogue of the Schwarzschild line 
element jTT] is 

ds'^ = -a^{z)df + dx"^ + dy^ + dz'^, (2.1) 

where z, x and y are analogues of radial r, (j) and 9 directions respectively. The 
lapse function vanishes at the event horizon placed at z = 0, and increases 
monotonically to unity as z increases from to cxd. 

2.1 Background Assumptions 

We consider two backgrounds with the specifications given below: 

1. Non- Rotating Background 

In this background, the fluid flow as well as the external magnetic field 
lines are along ^-direction (analogous to the radial direction) moving 
towards the black hole event horizon. The fluid velocity and the mag- 
netic field lines are taken to be V = u{z)e2 and B = Be^, respectively. 

2. Rotating Background 

This background demands the fluid to move in more than one dimen- 
sion. Thus we assume that the fluid motion and the external mag- 
netic field according to FIDO as V = V{z)e:x_ + u{z)ez and B = 
X{z)Bex + Be^ respectively, where B is a constant. 
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In MHD theory, plasma is treated to be a perfect fluid. We assume cold 
plasma surrounding the Schwarzschild black hole with equation of state 

— — — constant (2.2) 
Po 

which has vanishing thermal pressure and thermal energy. Here p is the 
mass density of the fluid. The fluid flow is assumed to be perturbed by the 
gravitational field of the black hole. The linearized perturbed quantities take 
the form 

p = pO + pp, V = V° + v, B = B° + Bb. (2.3) 

The variables with superscript zero are unperturbed quantities. The dimen- 
sionless notations for the perturbed quantities are 

p^p{t,z), v^v^{t,z)e^, h^b^{t,z)e^ (2.4) 

for non-rotating background whereas for rotating background 

V ^ v^{t,z)e^ + v^{t,z)e^, h ^b^{t,z)e^ + b^{t,z)e^. (2.5) 

We assume that the perturbed quantities have harmonic space and time 
dependence and hence these quantities can be expressed as follows 

Pit, z) = cie-^^-*-'^^), v^it, z) = C2e-^('^*-'=^), v,{t, z) = cae-'^'^*-'^^), 
h^{t, z) = C4e-'(^*-^^\ h,{t, z) = C5e-'('^*-'=^) (2.6) 

where Ci, C2, C3, C4 and C5 are arbitrary constants, k is the wave number and 
(jj is the angular frequency of waves. 



4 



3 Non-Rotating Background 



For the fluid flow in only one dimension, i.e., 2;-direction, we shall use the 
same Fourier analysed equations given in [20| as follows 

--C5 = 0, (3.1) 

a 

Lkc5 = 0, (3.2) 
Ci + Lku^ + C2 {(1 + -f^u^)Lk - (1 - 2-f\^) 

x(l + 7V)---7'^j =0, (3.3) 



u a 

ci7^{a^ + uu'C^ + l^u^)} + C2[7^(l + l^u"^) (^-^ + ^uk^ + 

WW{1 + l^u^){l + 47'm') + 2VaJ] = 0. (3.4) 

Equations (13. ip and (13.21) imply that C5 = which means that no perturba- 
tions occur in magnetic field of the fluid. Thus we are left with Eqs. (13. 3 p and 
(13.41) . It is mentioned here that for non- magnetized plasma, we also obtain 
the same two equations. 



3.1 Numerical Solutions 

We consider cold plasma of constant density. The time lapse is assumed to 
be tanh(2;) which vanishes at the horizon z = and increases to 1 as z — > cxo. 
Using mass conservation law in three dimensions, i.e., pju = constant we 
obtain the value of m = , ^ The determinant of the coefficients of 

Y'l+tanh^(2) 

the constants ci and C2 in Eqs. (13.31) and (13. 4p is a complex number. Solving 
this determinant by using the software Mathematical we obtain two complex 
values of the wave number k. Real and imaginary parts of this wave number 
show the propagation factor and the attenuation factor respectively for plane 
wave in dispersive medium. The propagation factor gives the phase velocity 
{vp = -j^^) and group velocity {vg = -jf^ij^) of the waves. The sinusoidal 
expressions then take the form g-'-M-^^i^-'-'^az) _ ^-L{ut-kxz)-k2z ^ where ki = 
Re{k) and k2 = Im{k). The wave numbers obtained are shown in Figures 1 
and 2. 

In Figure 1, the waves vanish in the region < 2; < 1.1 x 10~^ due to 
infinite wave number. The propagation factor of the waves is positive for 
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Surface Plot of Propagation Factor 



Surface Plot of Attenuation Factor 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 1: Wave propagation increases near the event horizon. Dispersion is 
found to be anomalous. 
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Surface Plot of Propagation Factor Surface Plot of Attenuation Factor 




Figure 2: The waves propagate rapidly near the event horizon. The whole 
region admits normal dispersion. 



the region < 2; < 1.1 x 10~*, < a; < 10. It increases with the increase 
in angular frequency and decreases with an increase in z. Thus the waves 
propagate rapidly far from the event horizon. The waves with higher angular 
frequencies move faster than the waves with lower angular frequencies. The 
attenuation factor increases with an increase in z. The waves damp as they 
move away from the event horizon. The phase velocity is less than the group 
velocity in the region < 2; < 10 which implies anomalous dispersion of the 
waves in this region |23j . There is an exceptional case of the waves with very 
small angular frequency that do not show this behavior. 

Figure 2 shows that the wave number is infinite in the region < 2; < 
1.1 X lO"*^. The propagation factor takes negative values in 1.1 x 10~* < 
2;<10, 0<ct;<lx 10~^^ which shows the same behavior as given in 
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Figure 1 except for this small region. The attenuation factor is negative in 
the region 1.1 x 10^^ ^ z < 0.68 and takes larger values aX z = 1.4 where 
the waves damp. The waves grow as the angular frequency decreases. The 
group velocity is less than the phase velocity for the region < a; < 10 which 
shows normal dispersion of waves. 



4 Rotating Non- Magnetized Background 

For the rotating background, the Fourier analysed Eqs.(4.13)-(4.15) of [20 
are 



CiL 



a 



+ ku] + C2 



LUJ 



a 
u' 



Yu + (1 + )ik 



-(1 - 2Yu^)n + Yu^)— + 2Yu^VV' 

u 



+ C37 



LUJ , 



a 



V + LkuV + u{{l + 2YV^)V' + 2YuVu'} =0, (4.1 



Ci72m{(1 + 7^2)^' + -^^uVu'} + C2i\ + ku^ fuV 

+{(1 + 27^2) (1 + 27^2) - -i^V^}V' + 27^(1 + 2-i^u^)uVu'] 

+C3[(-— + iku) 7^(1 + 7V2) + 7^u{(l + AfV'^)uu' 

+4T/y'(l + 7V^)}] = 0, (4.2) 
ci7^{a, + (1 + -i'^u^)uu' + 7^mVV'} 

+C2[7'(1 + l^u^)i (-^ + uk^ + 7'{«'(1 + 7'^i')(l + ^l^u^) 

+2uY{a, + (1 + 2Yu^)VV')]] + C37^ (-^ + uk^ uV 

WV'il + 47V2) +2V{a, + uu\l + 2Yu^)}\ = 0. (4.3) 



4.1 Numerical Solutions 

We consider the previous case with the modifications that V = tanh(2;), u = 
\ / T-n^^^TTTT ■ Solving the determinant of the coefficients of the constants 

y l+tanh^(z) ^ 

Ci, C2 and C3 in Eqs. fl4.ip - fl4.3l) . we obtain three complex wave numbers given 
in Figures 3, 4 and 5. 

In Figure 3, the wave number is infinite at the event horizon. The propa- 
gation factor decreases from the event horizon to ^ = 0.4 and then increases 
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Surface Plot of Propagation Factor 



Surface Plot of Attenuation Factor 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 3: Dispersion is normal for waves with negligible angular frequencies. 
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Surface Plot of Propagation Factor Surface Plot of Attenuation Factor 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 4: Many small regions admit normal and anomalous dispersion. 



to 2; = 6.25. It takes random values in the region 6.25 < z < 10. The atten- 
uation factor admits negative values in the region. The waves damp as they 
move away from the event horizon towards z = 0.55 and grow afterwards 
up to 2; = 4.5. The attenuation factor assumes random values in the region 
4.5 < 2; < 10. The phase velocity of waves is less than the group velocity for 
the region 0<2;<5, 0.1<a;<5 which implies anomalous dispersion. In 
the region < a; < 0.1, the phase is moving faster than the group of wavelets 
and the dispersion is normal. 

Figure 4 shows that no wave lies on the event horizon due to infinite wave 
number. The waves propagate rapidly with the increase in z and uj in the 
region 6 < 2; < 10. The attenuation factor increases with the increase in z 
and CO in the region < z < 4.25 which shows that the waves damp with the 
increase in z and uj. The waves with negligible angular frequency do not show 
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Surface Plot of Propagation Factor Surface Plot of Attenuation Factor 




Figure 5: The region admits normal dispersion for most of the points. The 
waves with neghgible angular frequencies show anomalous dispersion. 



this behavior. The phase velocity exceeds the group velocity for the regions 
0.55 < z < 3, 5<a'<10 and 1 < -2 < 3, 0.15 < a; < 5 which indicates 
normal dispersion. The group velocity is greater than the phase velocity for 
the region < z < 0.55, 5 < w < 10 where dispersion is anomalous. In the 
regions < z < 1, 0.15 < oj < 5 and < z < 0.15, there lie random points 
of normal as well as anomalous dispersion. 

In Figure 5, the wave number is infinite at the event horizon. The propa- 
gation factor increases with the increase in uj. It decreases in < z < 1 and 
increases in 1 < 2 < 6.3 with the increase in z and takes random values in 
the region 6.3 < 2; < 10. In the region < z < 4.6, the attenuation factor in- 
creases with the increase in u and z. Thus the waves damp with the increase 
in the values of angular frequency and z. For 0<2;<4.8, 5<ci;<10 and 
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1 < z < 4.8, 1 < < 2, the phase velocity is greater than the group velocity 
and the dispersion is normal. Rest of the region possesses some points of 
anomalous dispersion but mostly admitting normal dispersion. The waves 
with angular frequency between and 0.1 admit anomalous dispersion. The 
region 4.8 < 2; < 10 admits normal as well as anomalous points of dispersion. 



5 Rotating Magnetized Background 



For harmonic dependence of waves given by Eg. (12.61) . the Fourier analysed 
GRMHD Eqs.(5.15)-(5.20) of [20] can be written as 

€3(0;' + ika) — C2{{aXy + ikaX} + C5{aVy 

—C4{{auy — Luj + ikau} = 0, (5.1) 

C5 (^-^ + Lkuj = 0, (5.2) 

c^ik = 0, (5.3) 

c, f-^ + ,ku) + C2 l--7\ + ik{l + ^\^) + 2^\^VV' 
\ a / la 

-(1 - 27V)(1 + 7V)- j + C37' [(-- + Lku) V 

u ) l\ a / 

+u{{l + 27V^)y' + 2-f\Vu}] = 0, (5.4) 

cip^\{(l + 7V)F' + ^\Vu'} - —cMl - u^)ik + -(1 - u^) - uu'} 

An a 



+C2 



— I p7 uV — 



+ iku I p7 uV + 



An J a \ An J Ana 

+pf{{l + 27^2) (1 + 27V2) - 7V2}y' + 2p7^(l + 2-i'^u^)uVu'] 

+C3 - jpf (1 + 7V^) + ^y^ + cku jpf (1 + fy^) - ^ 



B'^ua' 
Ana 



+pru{{l + AYV^)uu' + 4(1 + YV^)VV'} 



0, 



(5.5) 



CiP7^[a^ + m{(1 + 7^M^)^x' + 7Vuy'}] 
+C2 - (p7^(l + 7^w^) + ^7^1 — + <! P7^(l + 7^M^' 



An 



a 



47r 



iku 



12 



+ {pj^u'il + 7'ti')(l + 47^1^2) + 2u^\{l + 2j\^)VV' + a,)} 



+— C4 X{l~u^)tk + X—{l-u^) - Xuu' +- — - =0. 




iku 



(5.6) 




Equations (15.21) and (15. 3p imply that cs is zero which gives that bz is zero. 
5.1 Numerical Solutions 

We have assumed the same velocity which we have taken in the previous 
section. In addition to this, we consider the same magnetic field components 
which we have used in the case of isothermal fiuid in [22], i-e., -B^ = BX = 



Under these conditions, when we solve the determinant of Eqs. ( 15TT]) . (15. 4p - 
(15.61) by using C5 = 0, we obtain a dispersion relation which leads to four 
complex wave numbers k. These wave numbers are shown in Figures 6-9. 

Figure 6 indicates that the wave number is infinite at the event horizon. 
The propagation factor has small values near the event horizon. It suddenly 
increases and then decreases as the waves move towards z = 1 after which it 
increases in the region < z < 6.3. This shows that the waves abruptly move 
faster, then their propagation decreases and afterwards increases as they 
move away from the event horizon in the mentioned region. The propagation 
of the waves is fast as we increase their angular frequency which is quite usual. 
The propagation factor takes random values in the region 6.3 < 2; < 10. 
The attenuation factor increases when the waves move away from the event 
horizon. It decreases in 0.3 < z < 0.7 and then increases monotonically in the 
region < 2; < 4, 1 < cj < 10. Thus damping occurs when the waves move 
away from the event horizon. The waves grow in the region 0.3 < z < 0.7 
after which they show smooth damping. 

The phase velocity is less than the group velocity for the waves near 
the event horizon which shows normal dispersion. As the waves move away 
from the event horizon, the group velocity overcomes the phase velocity and 
anomalous dispersion holds. The waves with low angular frequency admit 
anomalous dispersion. As the angular frequency increases, the waves show 
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Surface Plot of Propagation Factor Surface Plot of Attenuation Factor 




Surface Plot of Phase Velocity Surface Plot of Group Velocity 




Figure 6: Wave propagation is less near the event horizon. Anomalous dis- 
persion occur in most of the region. 
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Figure 7: Dispersion is random in most of the region. Negative wave propa- 
gation region is observable. 

normal dispersion near the event horizon but as these waves move away, they 
disperse anomalously. Thus these waves cannot transfer any information out 
of the region. 

Figure 7 demonstrates that the wave number is infinite at z = and 
thus no wave exists there. The propagation factor decreases abruptly when 
the waves move away from the event horizon. This factor increases with an 
increase in the value of z and uj. This means that the wave propagation 
abruptly decreases when the waves move away from the neighborhood of 
the event horizon of the Schwarzschild black hole. Then their propagation 
increases gradually. The waves with higher angular frequency propagate more 
rapidly in the environment than the waves with lower angular frequency. The 
attenuation factor decreases when moving away from the event horizon till 
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Surface Plot of Propagat ion Factor 



Surface P lot of Attenuation Factor 




Figure 8: Small regions of normal as well as anomalous dispersion are found. 
Negative wave propagation region occur. 

^ = 0.5. It then increases a little and decreases gradually in the region 
< z < 2, 2<a'<10. The waves grow as they move away from the event 
horizon, a little damping occurs after which they grow again in this region. 
In the rest of the region, the growth and damping occur randomly. 

The group velocity is greater than the phase velocity in the region 2 < 
2; < 10, 4 < a; < 10 which indicates anomalous dispersion. The rest of the 
region shows random points of dispersion. In the region 8 < 2; < 10, the 
phase as well as the group velocity admit the negative values. 

Figure 8 indicates that when the waves move away from the event hori- 
zon, the propagation factor decreases towards z = 1 and increases gradually 
towards z = 6. In the region 6 < 2; < 10, the propagation factor attains 
random values and in the region 0<2;<2, 2<a;<10, the attenuation 



16 



Surface Plot of Propagation Factor Surface Plot of Attenuation Factor 




Figure 9: Most of the regions contains random points of dispersion. Region 
with negative phase velocity propagation is also present. 

factor decreases as z increases. This indicates that the waves damp with the 
increase in z in this small region. In the rest of the region, the attenuation 
factor admits random values and thus the growth and damping of waves 
occur at random. 

The phase velocity is greater than the group velocity for most of the points 
in the regions 2 < 2; < 3, 0.5 < oj < 4 and 3<2;<4, l<a'<6 which 
indicates normal dispersion. The group velocity exceeds the phase velocity 
and the dispersion is normal in the regions (i) 1 < ^ < 2, 1 < uj < 10, (ii) 
2 < 2 < 3, 4 < w < 10 and (iii) < z < 3.5, 6 < < 10. The rest of 
the region contains points of normal as well as anomalous dispersion. The 
region 6 < 2 < 10 takes some points at which the phase as well as the group 
velocity have negative values. 
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In Figure 9, the wave number is infinite at ^ = 0. The propagation 
factor decreases at once on going away from the event horizon. It increases 
gradually when the value of z as well as the angular frequency increase in the 
region. This means that the propagation of waves decreases abruptly when 
they move away from the event horizon, and increases when the waves move 
faster. The waves with higher angular frequency propagate faster than the 
waves with lower angular frequency. The propagation factor assumes random 
values in the region 6 < z < 10. The attenuation factor has less variations 
in the region < 2; < 4 and shows increasing behavior as the value of z and 
ijj increases. This means that the waves damp either when they move away 
from the event horizon or when their angular frequency is increased. The 
variation of this factor is much more in the region 4 < z < 10. 

In the regions (i) < 2; < 3, 4 < a; < 10, (ii) l<z<4, 2<u;<4 and 
(iii) 0<2<2.5, 2<ij;<4, the phase velocity is greater than the group 
velocity and the dispersion is normal. The regions 3 < 2; < 4.4, 4 < a; < 10 
and 2.5<2;<5, 2<a;<4 possess anomalous dispersion due to the 
fact that the group velocity is greater than the phase velocity. The regions 
< ijj < 2 and 4.4 < z < 10 admit some points with normal and some with 
anomalous dispersion. In the region 5 < z < 10, the phase velocity takes 
negative values at some points and hence the region contains points with 
negative phase velocity propagation. 

6 Outlook 

In this paper, we work out the complex wave numbers near the event horizon 
of the Schwarzschild black hole. For this purpose, we have taken the pla- 
nar analogue of the Schwarzschild spacetime with non-rotating and rotating, 
non-magnetized and magnetized backgrounds. We have investigated |20] real 
dispersion relations for the same backgrounds by using Rindler coordinates. 
The Rindler coordinates give best approximation of the Schwarzschild space- 
time near the event horizon whereas this planar analogue is valid for the 
whole accretion cloud around the Schwarzschild event horizon. The proper- 
ties of the surrounding plasma are based on the complex wave numbers. 
The summary of the results obtained is outlined below: 

• The wave number is infinite at the horizon which supports the well- 
known fact that no information can be extracted from a black hole. 
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• There is no difference in the dispersion relations whether the plasma is 
magnetized or not, in the case of non- rotating background because the 
gravity produces no perturbations in the magnetic field. The fluid is 
dispersed normally in only one case (Figure 2) which allows the waves 
to pass through this region. The waves show anomalous dispersion in 
the other case. 

• The rotating background exhibits negative phase velocity propagation 
regions far away from the black hole event horizon. It seems that 
the black hole's gravity ceases the waves to admit the negative phase 
velocity propagation region. For the region with larger values of z, the 
dispersion is found to be normal at some points and anomalous at some 
other points. This creates a doubt whether all the waves can get out 
of this region. 

• There is only one case (Figure 5) in which normal dispersion occurs at 
most of the points in rotating non-magnetized background. 

• For the rotating magnetized background, many small regions which 
admit normal dispersion are found. These regions are surrounded by 
the regions of anomalous dispersion and hence no information can be 
extracted from these regions. 

The comparison of the results with those obtained from the previous 
literature is given as follows: 

• When we compare our results to that of [20|, we find that in non- 
rotating background we have detected the anomaly of waves which was 
not recognized previously. In the case of rotating magnetized back- 
ground we obtain the negative phase velocity propagation region far 
away from the event horizon whereas in [20], the whole region admits 
this property. 

• In the case of non- rotating isothermal plasma [21j, the dispersion is 
anomalous in all the cases whereas in cold plasma the region with 
normal dispersion is also found. It can be deduced that the plasma 
pressure ceases the chances of waves to come out of the region in the 
neighborhood of the event horizon in pure Schwarzschild background. 

• In the case of rotating magnetized isothermal plasma [22] , the medium 
admits normal dispersion of waves in plasma surrounding the event 
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horizon for the waves admitting high angular frequency. For the same 
case of cold plasma, the regions with normal dispersion are enclosed by 
the regions admitting anomalous dispersion of the waves. For the re- 
stricted Kerr black hole, it can be observed that the pressure reduction 
ceases the chances of waves to pass through the region. 

We conclude that no information can be extracted from the Schwarzschild 
event horizon. There are chances to collect the information from the exterior 
of the event horizon of the Schwarzschild black hole for the case of cold 
plasma. In the case of the restricted Kerr geometry, the cold plasma does 
not allow the waves to pass through in most of the cases. However, there are 
little chances for the waves to move out of the accretion disk. 
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